Modelling of the moving deformed triple 
contact line: influence of the fluid inertia 



Vadim S. Nikolayev 

ESEME, Service des Basses Temperatures, CE A- Grenoble, France and 
CEA-ESEME, ESPCI-PMMH, 10 rue Vauquelin, 75231 Paris Cedex 5, France 

Sergey L. Gavrilyuk 

University of Aix-Marseille & Polytech Marseille UMR CNRS 6595 IUSTI, 
5 rue E. Fermi, 13453 Marseille Cedex 13, France 

Henri Gouin 

University of Aix-Marseille & U.M.R. C.N.R.S. 6181, 
Case 322, Av. Escadrille Normandie-Niemen, 13397 Marseille Cedex 20 France 



Abstract 

For partial wetting, motion of the triple liquid-gas-solid contact line is influenced 
by heterogeneities of the solid surface. This influence can be strong in the case of 
inertial (e.g. oscillation) flows where the line can be pinned or move intermittently. A 
model that takes into account both surface defects and fluid inertia is proposed. The 
viscous dissipation in the bulk of the fluid is assumed to be negligible as compared 
to the dissipation in the vicinity of the contact line. The equations of motion and the 
boundary condition at the contact line are derived from Hamilton's principle. The 
rapid capillary rise along a vertical inhomogeneous wall is treated as an example. 
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1 Introduction 



In equilibrium, the wetting properties of a liquid in contact with a solid are 
well defined by the static contact angle 9 eq [1]. In the case of complete wetting 
(6 eq = 0), a thin prewetting film exists in front of the bulk of the liquid 
and its dynamics is well understood too. For partial wetting (9 eq ^ 0), a 
difficulty arises when the triple solid-liquid-gas contact line is moving. Huh 
and Scriven [2] showed that the viscous flow in contact line vicinity cannot 
be described by the no-slip condition, i.e. zero fluid velocity at the solid wall. 
Indeed, since the contact line belongs at the same time to the liquid and to 
the solid, the contact line velocity is ambiguous. This leads to a nonphysical 
divergence of the hydrodynamic pressure and of the viscous dissipation at the 
contact line. In reality, the dissipation is anomalously large but finite. 

A motion of the straight contact line is studied in the large majority of 
works. However, wetting dynamics experiments at partial wetting are almost 
inevitably submitted to an influence of surface heterogeneities that we call 
"defects" for the sake of brevity. On small scale, they lead to the contact line 
deformation; on macroscopic scale, to the contact angle hysteresis, i.e. to the 
difference between the advancing and receding contact angles. The contact 
line elasticity approach was proposed to describe the static deformation of the 
contact line by a localized defect [1]. It resulted in the logarithmic contact line 
shape, which was shown to describe correctly experimental data in the inter- 
mediate range of distances from the defect center. Obviously, the logarithmic 
shape fails to describe the contact line both very close to the defect and far 
from it where the contact line deformation needs to be finite. To obtain the 
correct description in the whole range, the influence of gravity [3] or the fluid 
mass conservation [4] needs to be taken into account. 

While the static hysteresis is understood relatively well [1,5,6], the contact 
line dynamics in the presence of defects is under active discussion, see [7] 
for related references. The contact line speed is often presented as a function 
of the dynamic contact angle 6 and of the static contact angle. However, an 
ambiguity appears in definition of this static value which should lie somewhere 
between the advancing and receding values of the static contact angle. 

Among theoretical studies of the dynamics in the presence of defects (that we 
call dynamics of the "deformed contact line" ) one can distinguish the papers 
[8,9]. Joanny & Robbins [8] proposed a general formalism for an arbitrary 
contact line shape but solved it only for the straight contact line. Golestanian 
and Raphael [9] used the quasi-static dynamics in conjunction with the contact 
line elasticity model to study the random deformation of the contact line. 

To describe the effect of defects, a single phenomenological parameter £, the 
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dissipation coefficient, was introduced [3,10] following [11]. £ contains all the 
physics of the contact line motion mechanism and has to be specific to the 
fluid-substrate system. To introduce it into the theory, one needs to assume 
that (i) the viscous dissipation in the bulk of the fluid is neglected with respect 
to that in the contact line vicinity, and (ii) the contact line dissipation is in- 
dependent of the direction of the contact line motion (advancing or receding). 
The energy dissipated per unit time at the contact line is then given in the 
lowest order of the contact line velocity v n (measured in the direction normal 
to the contact line) by the expression 

T = / ^ d/, (1) 



where the integration is performed along the contact line. The defects are 
modelled by the spatial variation of the surface energy of the support which, 
according to the Young-Laplace formula, is equivalent to the spatial variation 
of 9 eq . It was shown in [4] that in quasistatic approximation (i.e. when the 
fluid surface is assumed to be at equilibrium shape at any time moment) the 
expression ([I]) leads to the equation valid for any contact line (or fluid surface) 
geometry, 

v n = j (cos 9 eq - cos 9). (2) 

where a is the surface tension. Therefore, in any contact line motion model that 
results in (|2J), the coefficient £ can be interpreted as the contact line dissipation 
per its unit length. For instance, £ can be directly assimilated to the friction 
coefficient of Blake and Haynes [12]. The dissipation coefficient theory allows 
the contact line dynamics and the dynamic advancing and receding angles 
to be predicted once the surface heterogeneity and the initial conditions are 
defined [13]. 

The relationship ([2]) has been verified against multiple experiments (see [14] 
for a review). Using (jSJ), the £ value can be obtained from them. According to 
various experimental data, £ is much larger than the shear viscosity rj, which 
corroborates the assumption (i) above. However, there is still no certainty even 
about the order of value of ^/r? for a particular liquid-solid system. The exper- 
imental values range from 300 [11] to 10 7 [15] and depend on the conditions 
of experiment, the state of the substrate, etc. 

It is argued recently [7] that v n is non-linear in cos 9 which is equivalent to the 
statement that £ varies strongly with v n . At large v n , the observed deviations 
from the behavior (j2j) can be explained by the model [14] . However, one needs 
to treat the experimental £ value at small contact speeds with caution because 
it is usually determined [7, 15-17] from the slope of the dependence of (v n ) on 
(cos 9) (where the angle brackets mean averaging along the contact line) . It is 
shown in [13] that due to the collective effect of the defects, the dependence <^ 
with constant £ can result in a highly non-linear dependence of (v n ) on (cos 9) 
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near the pinning threshold where (v n ) is small. The inferred from the slope £ 
value turns out to be larger than its actual value. The influence of defects can 
thus lead to averaged v n {6) equations different from ([2]) at small contact line 
speeds. To separate the effect of defects from the physics mechanism of the 
contact line motion, a description of the dynamics of the contact line deformed 
by the defects is necessary. 

Most theoretical approaches aim to describe the slow fluid motion and thus 
neglect the fluid inertia. However, the inertial effects can be important in many 
common situations. One can mention the early stages of drop spreading [18], 
possibly after its impact with a support [19], the capillary rise [20-22], and 
oscillating flows in various geometries: those of sessile or pendant drops [23-27], 
liquid in containers [28-31], or capillary bridges [32]. 

In the partial wetting regime all these situations involve a period of time 
where v n is small and then the defect impact on dynamics is large. Indeed, v n 
necessarily passes through zero during oscillations. As for relaxation processes, 
the influence of defects becomes important at the end of the evolution. 

This problem is especially important for oscillation flows where the choice 
of the boundary conditions at the solid walls is a longstanding problem [28]. 
Evidently, when the amplitude of oscillations is small enough, the contact 
line is pinned [24], so that the fixed position boundary condition need to be 
applied. When the amplitude is larger, the contact line can be pinned during 
a part of the oscillation period [29,30]. At large amplitude, it moves almost 
all the time. In the relatively well studied [26, 31] first regime, there is no 
singularity at the contact line. The transition from the first to the second 
regime and especially the complicated dynamics in the second regime did not 
yet find its explanation. In this paper we develop a model suitable to solve 
these problems. 

The equations of motion are derived in the "shallow water" approximation in 
sections [2] and [3J This approximation permits to solve analytically a model 
problem (section S]) in order to validate the model and find a first-order "in- 
ertial" correction to the quasi-static behavior. The results are summarized in 
the section [5j 



2 Equations of motion 

The oscillating flows are considered often in the inviscid fluid approxima- 
tion [27,31]. This is justified by the thinness of the viscous boundary model at 
moderate Reynolds numbers where the flow is not yet turbulent. For the con- 
tact line problems, an additional argument related to the anomalously strong 
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Fig. 1. Sketch of the liquid layer near the vertical (possibly moving with the speed 
v) wall with surface defects. 

contact line dissipation can be put forward. Indeed the experiments [24,28-30] 
show that its contribution can be larger than that of the viscosity in the liquid 
bulk. Therefore, we will neglect the viscous dissipation in the bulk by assum- 
ing that the liquid is inviscid. The only source of dissipation is assumed to be 
at the contact line. 

A Cartesian reference frame Oxyz is used to describe the fluid motion in the 
domain x > 0, — oo < y < oo, < z < h(t,x,y) (Fig. [1]). The free surface is 
described by the equation z = h(t, x, y) and is assumed to be weakly deformed 



where the gradient is taken with respect to the variables {x,y). The vertical 
wall (x = 0) may move in z direction. The gravity is directed downward. The 
contact line is described by the equation 



For convenience, the function h(t, x, y) is assumed to be periodic (with the 
period 2L) in the ^/-direction. This is not a restrictive assumption since one 
might put L — > oo in the resulting equations. 

To derive the governing equations and boundary conditions in the absence of 
the volume dissipation, we use the Hamilton principle of stationary action. 
The idea is to formulate the Lagrangian, and to present the variation of the 
Hamilton action as a linear functional of virtual displacements allowing to find 
not only the governing equations but also the natural boundary conditions 



In the present work, we will use the "shallow water" approximation which will 
allow the analytical results to be obtained and analyzed. The Hamilton action 
a and the Lagrangian £ (which is the difference between the kinetic and the 
potential energy of the system) read 



\Vh\ < 1, 



(3) 



z = h w (t,y) = h(t,0,y). 



(4) 



[33,34]. 
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a r°° f L { i \ a r L rh w (t,y) 

TlL dl / i (v /rT W- 1 ) d! ' + 2l/- l d U c( «- z)iz - (5) 

Here p is the density of the fluid assumed incompressible, u = (u x ,u y ) is the 
fluid velocity field. The first integral in the Lagrangian is a difference between 
the kinetic and the gravitational potential energy of the system written in the 
shallow water approximation. The second integral is the liquid-gas interface 
energy, and the last integral is the energy of the fluid-solid surface [5], crc(y, z) 
being the difference between the interfacial tension of gas-solid and liquid-solid 
interfaces. According to the Young expression, c(y,z) = cos(9 eq ), where 9 eq is 
the value of the equilibrium contact angle. The variation of the latter along 
the substrate (i. e. along the y — z plane) reflects the surface heterogeneity 
(surface defects). 

The variation of the action is taken under the mass conservation constraint 

^ + V-(hu) = (6) 
and the non-penetration boundary condition at the vertical wall 

u*\ x =o = 0- ( 7 ) 



Search of the extremum for the action a (see Appendix |A} results in the 
following fluid motion equation: 



d (phu) 
dt 



+ V • (phu g) u) + pghVh - ahV(Ah) = 0. 



Its back substitution into the Lagrangian ([5]) reduces the latter to the form 
(see also flA~3l) ) 



5C 



a 





'dh 


L 


dx 



+ c(y, K 



2=0 



5h w dy. 



(9) 



The dissipation is introduced [3] in the usual way through the dissipation 
functional, (CQ). For this, we need to find an expression for the normal velocity 
v n of the contact line. When the vertical wall moves vertically with respect to 
the reference system with the velocity v, the velocity of the contact line with 
respect to the wall reads 



' dh u 
~dt 



+ v 



dh w 



-1/2 



dh u 
~dt 



+ v, 
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where v is assumed to be positive for the downward plate motion. Since the 
volume dissipation inside the fluid is neglected, the dissipation functional does 
not depend on the bulk fluid velocity u and the dissipation account does not 
influence the governing equation (|SJ). 

Following [3], one can write now the generalized Euler-Lagrange equation 
which takes into account the dissipation term 



d ( 5C 



dt V 5h, 



5C 

Sh„, 



ST 
6h,„ 



(10) 



where 5 . . . /5 . . . means functional derivative and h w denotes the time deriva- 
tive. The substitution of ([9]) into ( fTOi) results in the equation 



a 



dh 



dx 



+ c(y, h w ) 



x=0 



dh w 
dt 



that defines the contact line position h w . Eq. (fill plays a role of the boundary 
condition for the momentum equation (JHJ). 

Let us note that the shallow water approximation is not necessary to obtain 
the main result (1111) . It is known [34] that if the full Lagrangian (involving 
the integration over z of the kinetic energy) is used instead of its "shallow 
water" counterpart (jSJ), the classical Euler equation of motion for the inviscid 
fluid result is obtained instead of the equation ([8]). The back substitution into 
the Lagrangian would lead to the same expression (Q, which is necessary to 
obtain Eq. (fTTD. 



3 Reduction of governing equations 



The boundary conditions for the governing equations (0 EJ) should be com- 
pleted by the condition at infinity 

/il^oo = d, (12) 

where d is the constant depth. The following initial conditions will be used, 

u\t=o = 0, h\ t =o = d. (13) 

The governing equations written for the deviation from the initial state h = 
h — d can be linearized: 

Oh du 

— + dV ■ u = 0, p— + pgVh - aVAh = 0. (14) 
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The velocity u is eliminated from these equations by applying the time deriva- 
tive to the first of them and subtracting the divergence of the second. The 
result reads 



P 



d 2 h 



- pgdAh + adA 2 h = 0. 



(15a) 



Evidently, the first time-dependent term of this equation is inertial. Two other 
terms correspond to the gravity and to the surface tension contributions re- 
spectively. Without the first term, Eq. (j!5aj) would reduce to the quasistatic 
linearized equation for the interface shape [3]. 



The boundary and initial conditions for (j!5aj> are obtained from (0 fT2HT3"j) by 
using (CO 



h\ 



(pgh — a Ah 



dx 



0. 



x=0 

h\ t =o = 0, 
dh 



dt 



0. 



(15b) 
(15c) 

(15d) 
(15e) 



The nonlinear equation (II ip rewritten for h 



dh 



dh 

dx 



+ c(y, h w ) 



x=0 



(15f) 



where h w = h w — d, closes the problem statement. 



Only the solution for h w (t,y) is of interest. However, it can only be obtained 
by analyzing the full problem (fl5l) for h. The set of equations (jT5l) allows h w 
to be obtained for arbitrary distribution of the surface defects c(y, z). 



4 Solution of a model problem 



To demonstrate the applicability of the developed formalism, let us solve a 
problem of the capillary rise on an immobile (v = 0) wall with a single stripe 
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defect (see Fig. [3] below) 



J c d , 


\y\ 


{ c„ 


\y\ 



c(y) = { ' - - (16) 



where q, c s and w are constants, i.e. for the simplest case in which c is inde- 
pendent of z. The values of q and c s are the cosines of the equilibrium contact 
angles inside and outside the stripe defect respectively, w being its half-width. 
This problem has already been solved in quasi-static approximation [3] (in 
which the fluid motion is neglected). We choose the same setup to show the 
role of the fluid inertia by comparing the quasi-static and inertial solutions. 

Since z dependence in (JTB1) is absent, the problem (fTS]) becomes linear. Fol- 
lowing [3], h and h w are broken into parts: 

h = d + h {t,x)+h 1 (t,x,y), 
h w = d + h w0 (t) + h wl (t,y), 

where h and h w0 are the averaged over the y-direction values of h and h w 
respectively. The problems for h$ and h\ are decoupled and can be solved 
separately. 



4-1 Solution for h 



wO 



Consider first the problem for ho that can be obtained from (fl5|) by accounting 
for the fact dho/dy = 0. By applying the Laplace transform 



h (x) = / ho(t, x) exp(— st) dt, (18) 
Jo 

one obtains from (\15h .d.e) the ordinary differential equation 

-^-20^ + ^ = 0, (19) 

where a = pg/(2a), [3 = ps 2 /(ad). The general solution of (|T9l) has a form 

h (x) = ai exp(-pix) + a 2 exp(-p 2 x) + a 3 exp(pia;) + a 4 exp(p 2 ^), (20) 
where 



Pi, 2 = da ± J a 2 - (3, Rep li2 >0. (21) 
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The boundary conditions for ho have the same form (115b . c) as for ho', (115f[) 
transforms into 



sh w o 



a ( dhn 



£ \ dx 



+ 



x=0 



where h w o is the Laplace transform of h w0 . Because of ( Il5bl) . 

a 3 — 04 = 0, 
and (j!5cp results in the following condition 
dh 



dx 



-h w0 (pi +P2) = -h w0 \ 2 [a + J 13 



x=0 



(22) 



(23) 



(24) 



The final form for h w o can be obtained by the substitution of f l24"l) into ( |22l : 

-1 



a 



hwo c s 



s s + r„ 1 \/sv / 2^ + 1 



(25) 



where r r = Z c £/cx is a characteristic time scale for the quasi-static relaxation [3] 



and Ti = l c j 1 \fgd is the characteristic inertial time, l c = Ja/(pg) being the 
capillary length. The inverse Laplace transform of ( l25i) can be found rigorously 
as described in Appendix. It reads 




2VT+7 2 



exp 



— (r + VTT 



exp 




— (r-VTT 



(26) 



where erfc(-) is the complementary error function [35], and r = Ti/r r char- 
acterizes a relative importance of the fluid inertia. It is related to the Weber 
number which can be obtained from fU5al) where both the terms of inertia 
and of surface tension are present. Taking r r and l c as the scales of time and 
length, one can obtain for these terms ph/r 2 and adh/lf. Their ratio is the 
Weber number, We = r 2 . 



In the limit r — >■ 0, 



reduces to the quasi-static result [3] 
h w o(t)\ r ->o = c s l c [l - exp(-t/r r )]. 



(27) 



The difference between these two dependencies is transparent from Fig. [2j the 
fluid inertia slightly slows down the contact line relaxation. For r < 0.01, the 
difference between the functions ( )26|) and (|27|) is very small so that in practice 
the exponential solution ( 1271) is valid. 
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1 2 3 4 5 6 

t/T 
r 

Fig. 2. The time dependence of the y-averaged hight of the contact line rise h W Q. The 
exponential quasi-static law (|27|) (solid line) is compared to the inertial behavior 
(f26l) (dotted line) calculated for r = 0.5. 

4-2 Solution for h w % 



The statement of the problem for h\ is given by (TT5l) with an additional condi- 
tion of zero y- averaged value of h\. For a single defect, this condition reduces 
to h\{y — > ±oo) =0, which allows the Fourier transform 



h x {t,x) 



h 1 {t,x,y)exp(-iky) dy 



(2* 



to be applied. The Laplace transform ffT8]) can further be applied to hi(t,x) 
and results in hi(x). 

The further procedure is fully analogous to that developed in the previous 
section. It leads to the following equation for h w i 



- o ( dhi 
sh w i = — — — 
£ \ ax 



x=0 



+ H1 

S . 



(29) 



where C\ is the Fourier transform of c(y) — c s , 

sin kw 
c\ = 2oc- 



k ' 

and 5c = Cd — c s . The hi derivative is related to h w i through 



dhi 



dx 



—h w \ 



PiP2(pi +P2) 



x=0 k 2 +p lP2 
where pi^ are defined by (|2T|) . in which 

P9 , ,2 n P 2 . P9,2 . >4 

a — - — h k , p = — -s H k + k 



2a 



ad 



a 



(30) 



(31) 



(32) 
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should now be used. The substitution of (!3T!) into (|29|) leads to an expression, 
for which the inverse Laplace transform is difficult to find. To overcome this 
difficulty, one can make use of the smallness of the parameter r, i.e. assume the 
smallness of the contribution of the inertial effects. In this limit, (l3"Tj) reduces 
to 

dhi 



dx 



~h w i 



2 = 



k 2 + lz 2 + B(rT r s) 



where the coefficient B is 



B 



k 2 + It 2 k 2 + 1.51: 



\k\ 



k 2 + I- 2 



(33) 



(34) 



Note that B is positive and diverges at small k like B ~ \k\ 1 . The equation 
for h w i becomes 



h w \ — c\ 



s Jk 2 + l~ 2 + l-\s + B(rT r s) 



n -1 



(35) 



This expression is rational in s and its inverse Laplace transform can be found 
using conventional methods: 



h w i(t) — C\ 
Here 



1 



1 



k 2 + L- 2 B(rr r ) 2 ( Sl - s 2 ) 



exp(sit) exp(s 2 ^) 



si 



«2 



(36) 



s l,2 



1 ± ^ll-ABl c r 2 Jk 2 l 2 + l 



2Bl c r 2 r r ^ 
To obtain the time-varying contact line shape h w i(t,y), one needs to apply 
the inverse Fourier transform to (I3"6"|) . It can be done numerically. The result 
is shown in Fig. [3J 

At r — >• 0, one recovers the quasi-static result [3] 



h w i(t)\ 



Cl 



k 2 + 1: 



1 - exp Jk 2 l 2 c + 1 



(3J 



which corresponds to a non-zero value at k — > 0. However, by considering the 
small /c limit of (1361) , one arrives to the result 



h w i(t)\ 



Cit 2 

2BV 2 



0. 



(39) 



This small k ambiguity leads to a small but regular numerical error in h w i(t, y) 
when the inverse Fourier transform of ( |36|) is found at very small r. In practice, 
the expression ( 1381) should be used instead for r < 0.05. 



1 Notice that although s\ t 2 become complex at small k, h w \ remains real. 
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Fig. 3. The contact line distortion h w \ calculated for 



0.5. The defect area 



(half-width of which is w = 0.2/ c ) with cosd eq = a is shadowed. The rest of the 
plane has cosO eq = c s . The parameter of the curves is the time t in the r r units. 
In the insert, the result of the inertial theory for r = 0.5 (solid curve) is compared 
to the quasi-static (r — ► 0) result (dotted curve). Both curves are calculated for 

t = 0.17V. 

Similarly to the case of h w o, the inertial slowing down manifests itself in the 
case of h w \. The insert in Fig.[3]shows that the topmost point of the contact line 
rises slower when the hydrodynamic motion is taken into account. However, 
the hydrodynamic effects are not limited to the simple slowing down, they 
change the shape of the contact line deformation. In the beginning of the 
capillary rise, the joint effect of the inertia and fluid mass conservation creates 
two wells in the h\ shape (Fig. EJ) from which the fluid flows off to form a 
bump in the middle. Later on, these cavities become more and more shallow 
and wide until they disappear in the equilibrium static contact line shape [3]. 



5 Results and discussion 



The main result of the present paper is a model that allowed to describe the 
inertial regimes of the contact line motion whenever the impact of surface 
defects is important in the partial wetting regime, e.g. the pinning-depinning 
of the contact line during its oscillation. The equations of motion were derived 
for the plane vertical wall geometry and can be applied for the arbitrary defect 
pattern and both for the spontaneous and forced contact line motion. The 
hydrodynamic shallow water approximation was used to obtain the solution 
in the closed analytical form. However, this approach can be easily generalized 
to any geometry and fluid depth. 

Indeed, let us have a closer look to the equation (fTTI) . Its most interest- 
ing feature is its invariance. In the hydrodynamic theory, it has exactly the 
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same form as in the quasi-static theory [13]. Since c(y,z) = cos[9 eq (y, z)\ and 
dh/dx\ x=0 ~ — cos 9 under the assumption d3J), one notices that (iTTj) is nothing 
else than the equation ([21 . It was derived under the assumption of the small 
deformation of the straight contact line. By analogy with the quasistatic re- 
sult [4], it is however reasonable to assume that (T5J) is valid for any contact line 
deformation. Any contact line inertial problem can thus be solved by applying 
the Euler set of equations for the inviscid fluid with the boundary condition 
(j2J) where 9 eg (y, z) models the surface heterogeneity. 

The asymptotics of the capillary rise at small times which follows from the 
present model can be compared to the available experimental data. It stems 
from the equation (J2J) that v n is constant at t — because 9(t = 0) is fixed by 
the initial condition for h. Since h(t = 0) = has been chosen, 9(t = 0) = tt/2 
and h(t) ~ t at small times independently of the geometry of the problem. This 
asymptotics agrees with the results [21] on the capillary rise between flat walls. 
The asymptotics h(t) ~ t 3//2 obtained for the complete wetting case [20, 36] 
cannot be compared to our model restricted to the partial wetting. As a matter 
of fact, due to existence of the prewetting film [1], the contact line dissipation 
anomaly is much weaker for the complete wetting. The kinetics of the capillary 
rise is thus defined by the balance of the surface tension and inertia [20] which 
results in h(t) ~ t 3 / 2 . 

The effect of inertia on the contact line relaxation is two- fold. First, it is the 
slowing down of the average contact line motion. When the contact line speed 
is inhomogeneous, the inertia amplifies its variation along the contact line. 

The impact of the inertia on the contact line relaxation can be measured by 
the value r = Ti/r r of the ratio of two characteristic times. The first of them is 
the inertial time and the second is the quasi-static relaxation time. The quasi- 
static result is recovered in the limit r — > 0. Note that the Weber number (the 
ratio of the inertial and surface tension terms) is related to r: We = r 2 . 

The simplicity of the developed approach is its main advantage. It uses a single 
phenomenological parameter (the contact line dissipation coefficient £) and 
allows complicated three dimensional problems to be treated. In this article 
we analyzed the contact line motion in presence of a stripe defect. This model 
can be directly verified against experiments [37-39]. However, the information 
presented in these papers is insufficient for a direct comparison and additional 
measurements are necessary. 
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A Derivation of the governing equation ([8]) 



Assumption ([3]) permits us to rewrite (jSJ) in a simpler form 



c = -Tdx [ L f^M 

2L 7o i-L V 2 



2L Jo 



L \Vh\ 



a 



2L J-l Jo 



2 



dy - 



c(y, z) dz 



(A.1) 



The variation of the Lagrangian (1A.1I) submitted to the constraint (jSJ) is taken 
in several steps. First, we introduce a smooth one-parameter family of virtual 
motions 

&=$(t,X,e), x = {x,y) 

(X stands for the Lagrangian coordinates, e is a small parameter in the vicinity 
of zero and x = <p(t,X) = <&(t,X,0) is the real motion). Then, we define the 
virtual displacements as functions of (t,X): 



5x(t, X) 



d$(t,X,e) 



de 



£=0 



Later, we consider the virtual displacements as functions of (t, x) due to the 
possibility to invert X through X = The variations of the free sur- 

face h(t, x) and the velocity field u(t, x) compatible with the mass conservation 
law (jSJ) are given in the motion space (t,x) by (see, for example, [34] or [40]) 



5h(t, x) = —V • (hSx), 8u(t, x) 



ddx du 



dt dx 



5x. 



(A.2) 



u ■ V 



Here 

d_ _ d_ 
dt~ dt 

is the material time derivative. Finally, the expression of the variation of a in 
terms of the virtual displacements reads 



2L5a 



t2 poo pL 

dt dx 

ti JO J-L 



o IxTI / cLSx du \ 

V-(kS^ + P ku.l—--J 1 U 



pghV ■ (hSx) 



rt2 roo plj 

dy + a dt dx Vh ■ V [V • (hSx)] dy + 

Jti Jo J-L 

rt 2 rL 

a dt c(y,h w )5h w dy. 
Jti J-l 
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By integrating by parts, accounting for the periodicity in y-direction, and the 
boundary condition ((Tj) in the form 



Sx ■ n\ 



x=0 



0. 



where n = (1, 0) is the unit normal vector to the plane x = 0, we obtain 



2L5a 



roc pL 

dt dx 

ti JO J-L 

a 



d (phu) 
dt~ 



V • (phu ®u)- hV(pgh) + ahV(Ah) 



5xdy + 



rl2 roc rt2 rlj 

/ dt / dx V- [V- (h8x)Vh)dy + a / dt / c(t,y,h w )5h w dy, 

Jti JO J-L Jti J-L 

where ® denotes the direct product of the vectors. The definition of 5h w 

8h w = - V • (h8x)\ x=0 
gives us the final expression for the variation of the Hamilton action: 



rt2 roo rL 

2L5a = dt dx 

Jt x Jo J-L 



d (phu) 



-ah\7(Ah) 



dt 

rt 2 



V • (phu <S> u ) — hV(pgh) 



rt2 rJ^ 

5xdy + a dt [Vh ■ n + c(y,h w )}5h w dy. (A. 3) 
Jt x J-L 



The variations 5x and 5h w are independent and the term containing 5x is 
absent in the dissipation function because of neglect of the viscosity in the 
bulk of the fluid. Therefore, the integral containing 5x in (IA.3I) has to be 
equal to zero, which results in the governing equation ([H]). 



B Laplace-inversion of ( 125 



The expression ( 1251) has a form 



f( s ) = [ s (s + 2bVs + m 2 )}-\ 
for which the inverse Laplace transform is sought. The inverse of 



(B.l) 



f(s — m 2 



where 



1 



2m?b 



1/2 



+ 



1/2 



U 2 {U2 - Uij 



yfs — m yfs + m yfs — Mi 



u x [ui - u 2 ) 
\fs~u 2 



(B.2) 



u 1>2 = -b± Vb 2 + m 2 , (B.3) 
is found using the inverse of (^/s + m) _1 listed in the tables [35]. To find the 
inverse of f(s), the original of f(s — m?) should be multiplied by exp(— m 2 t) 
according to the rules of the Laplace transform. 
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